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Interfacial superconductivity is observed in a variety of heterostructures composed of different materials 
including superconducting and nonsuperconducting (at appropriate doping and temperatures) cuprates 
and iron-based pnictides. The origin of this superconductivity remains in many cases unclear. Here, we 
propose a general mechanism of interfacial superconductivity for systems with competing order parame¬ 
ters. We assume that parameters characterizing the material allow formation of another order like charge- 
or spin-density wave competing and prevailing superconductivity in the bulk (hidden superconductivity). 
Diffusive electron scattering on the interface results in a suppression of this order and releasing the super¬ 
conductivity. Our theory is based on the use of Ginzburg-Landau equations applicable to a broad class 
of systems. We demonstrate that the local superconductivity appears in the vicinity of the interface and 
the spatial dependence of the superconducting order parameter A(x) is described by the Gross-Pitaevskii 
equation. Solving this equation we obtain quantized values of temperature and doping levels at which A(x) 
appears. Remarkably, the local superconductivity shows up even in the case when the rival order is only 
slightly suppressed and may arise also on the surface of the sample (surface superconductivity). 

PACS numbers: 74.70.Xa, 74.62.Dh, 74.20.De, 74.78.Fk 


I. INTRODUCTION 

Interesting phenomena have been discovered few years 
ago in the study of superconductivity in different ma¬ 
terials, especially in high-Tc superconductors—cuprates 
and Fe-based pnictides. It turned out that the crit¬ 
ical temperature of the superconducting transition 
in heterostructures, e.g., in bilayers, is higher than the 
critical temperature Tc of bare films that can be even 
nonsuperconducting. The authors of Ref. 14 (see also 
Ref. 9) studied bilayers consisting of two cuprates— 
overdoped (La2-j:Sr;cCu04 with x = 0.35) and underdoped 
(0.1 <x< 0.12) cuprate films. The largest increase of Tc 
was about 11 K—from Tc = 21 K in bare underdoped films 
to Tc = 32 K in bilayers. Moreover, superconductivity has 
been observed in bilayers composed of nonsuperconduct¬ 
ing materials, La2-j:SrjcCu04 with x = 0 (an insulator) and 
La2-jcSr;cCu04 with x = 0.45 (normal metal). ^ The critical 
temperature reached 50 K. Further experimental studies 
showed the independence of the critical transition temper¬ 
ature in bilayers La2-jcSrj:Cu04/La2SrCu04 on x in a rather 
wide doping interval (0.15 < x < 0.47).^^ 

A slight increase of Tc has been measured also in an¬ 
other high-Tc superconductor—YBa2Cu307 films covered 
by a thin Ag film.^® A few decades ago, a similar effect has 
been observed in bilayers composed of conventional low- Tc 
superconductors and normal metals.These latter ex¬ 
periments were motivated by original Ginzburg’s ideas on 
the possibility to get a surface superconductivity with a high 
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transition temperature. ’ 

Another high-Tc superconductors where an enhance¬ 
ment of superconductivity at the interface has been es¬ 
tablished are the so-called iron-based pnictides discovered 
in 2008.^^ In this type of superconductors competing order 
parameters (OP), magnetic and superconducting, may co¬ 
exist. To be more exact, the spin density wave (SOW) and 


the superconducting OP may arise in these materials (see 
reviews Refs. 23-27) with the amplitudes 14/ and A depend¬ 
ing on temperature and doping. The interfacial supercon¬ 
ductivity in one of the Fe-based pnictides (CaFeAs doped 
with La, Ce, Pr or Nd) was observed by Wei etal?^ Whereas 
the bulk critical temperature was equal to about 30 K, a 
small fraction of samples had a Tc - 49 K. Even higher 
Tc - 77 K was achieved by a Chinese experimental group 
in another Fe-based material (single unit-cell FeSe films on 
SrTiOs).^^ 

Very encouraging for understanding the very nature of 
high- Tc superconductivity seems to be the effect of appar¬ 
ent enhancement of superconducting transition tempera¬ 
ture at the interface between an iron-based chalcogenide 
superconductor (FeSe) and SrTiOa used as substrate.^®’^^ 
This discovery questioned the role of phonons in bulk 
iron-based superconductors^® and confirmed the pres¬ 
ence of the magnetic order (spin-density wave) as a key 
ingredient for high-Tc superconductivity in iron-based 
superconductors.®^ 

Rather actively the interface superconductivity is stud¬ 
ied in heterostructures LaAlOs/SrTiOs.® It is assumed that 
a two-dimensional electron gas is formed at the interface. 
Effect of an electric field has been employed to explore the 
phase diagram of LaAlOa /SrTiOs interface.®’"^ Aside from su¬ 
perconductivity, the phenomenon of ferromagnetism in¬ 
duced at the interface of an oxide heterostructure has been 
observed recently.®®’®® The review in Ref. 34 provides an 
excellent overview over the possible symmetries and de¬ 
grees of freedom of correlated electrons that can evolve at 
oxide interfaces. It includes, inter alia, superconductivity, 
magnetism, ferroelectricity, and charge- and spin orders as 
well. Moreover, it has been found, that at the interfaces be¬ 
tween LaAlOs and SrTiOs, superconductivity coexists with 
ferromagnetism,®®'®®"®® a surprising result offering a poten¬ 
tial for exotic superconducting phenomena due to highly 
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broken inversion symmetry of the interface and a ferromag¬ 
netic background.^® 

Several theories have been suggested to explain the phe¬ 
nomenon of interface superconductivity. Some of them 
consider a nonuniform charge distribution near the inter¬ 
face and use a phenomenological relation of Tc to this 
distribution.^’*^®’®® Different ideas have been used in other 
theories,®® where bilayers composed of superconductors 
with different ratio of the Tc and pairing strength were con¬ 
sidered and it was assumed that Tc is suppressed by phase 
fluctuations. In the vicinity of the interface the role of these 
fluctuations is not important as compared to the bulk due 
to suppression of fluctuations by the proximity effect. How¬ 
ever, this suggestion cannot explain the observed indepen¬ 
dence of Tc on the doping level x. Perhaps, there is no single 
mechanism responsible for the interface superconductivity 
because it was observed in quite different materials under 
various conditions. The important ingredient of this effect 
is the presence of an interface or some sort of nonhomo¬ 
geneity. 

In this paper we propose a new mechanism for the inter¬ 
face superconductivity. The proposed mechanism is very 
robust and general being independent of the microscopic 
details of considered materials. It is applicable to any ma¬ 
terials where, alongside with the superconducting order pa¬ 
rameter (OP), another OP exists. We do not pretend to apply 
our theory to any system where the interface superconduc¬ 
tivity occurs, but we show that it can be used to materials 
in which two OPs may potentially exist. As is well known, 
in high- Tc superconductors, an important role is played by 
the charge- or spin-ordering.®®”®'" In cuprates, a charge den¬ 
sity wave (CDW) has been observed recently in numerous 
experiments^®”®® and discussed in Refs. 51-55. It may ex¬ 
ist alongside with superconductivity, whereas in Fe-based 
superconductors, the spin density wave is more impor¬ 
tant. The presence of the nonsuperconducting OP W (CDW 
or SDW) changes such characteristics of superconducting 
state as London penetration depth,®®”®® heat capacity,®®’®® 
etc. 

We show that the presence of a ‘hidden’ OP A in a het¬ 
erostructure with W iiO may lead to appearance of local su¬ 
perconductivity at the interface at temperatures T > TdW), 
where TdW) is the critical temperature of the supercon¬ 
ducting transition of bare films composing the heterostruc¬ 
ture which depends on the amplitude W. We consider a 
heterostructure with an interface where W is locally sup¬ 
pressed. This suppression may be caused by an enhanced 
impurity scattering and doping level in the vicinity of the in¬ 
terface. The enhanced impurity scattering can be caused by 
the interdiffusion of atoms and/or roughness of the inter¬ 
face. Both factors suppress the CDW or SDW. In this case, 
in a vicinity of the interface where W is suppressed, lo¬ 
cal superconductivity arises with A(x) decaying on a char¬ 
acteristic scale of the order of the superconducting coher¬ 
ence length fs- Interestingly, the local superconductivity 
occurs at ‘quantized’ temperatures because the OP A(x) is 
described by the linearized Gross-Pitaevskii equation, i.e., 
by the Schrodinger equation with a one-dimensional po¬ 


tential well which always has discrete energy levels (or only 
one level). In one-dimensional case (flat interface) the local 
superconductivity arises even at a rather small suppression 
of W. 

Note that stimulation of the bulk superconductivity by 
impurities in materials with two OPs was considered by 
one of tbe authors a long time ago.®* Recently, the effect of 
superconductivity stimulation by impurities in tbe bulk of 
Fe-based pnictides with two OPs has been analysed by Fern- 
nades et al. in Ref. 62. However, to tbe best of our knowl¬ 
edge, there has been no study of the mechanism of the in¬ 
terface superconductivity presented in this paper. 

In Section II, the general nonhomogeneous Ginzburg- 
Landau equations describing two coupled order parameters 
are introduced. In tbe homogeneous case, we derive con¬ 
ditions for the coefficients in the Ginzburg-Landau equa¬ 
tions that should be fulfilled for obtaining one of the possi¬ 
ble three phases: 1) pure superconducting, 2) pure charge- 
or spin-density wave, 3) a mixed state. In the nonhomoge¬ 
neous case, a solution of the Ginzburg-Landau equations 
for the order parameters is obtained, whereby detailed cal¬ 
culations are shifted to the Appendix. In Section III, we dis¬ 
cuss applicability of the theory to the case when the con¬ 
stituents of the heterostructure are high- Tc cuprates or iron- 
based pnictides. We propose also an experimental setup 
suitable for testing our predictions. Concluding, we discuss 
our results in Section IV. 


II. GENERAL CONSIDERATIONS 


A. Free energy and self-consistency equations 


We start with the expression for the free energy ^ for a 
system with two order parameters (OPs), A and W. In one¬ 
dimensional case, i.e., in case of a preferred direction pro¬ 
vided by the interface, the free energy is given by 
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dx 


(AT - flsA^-t -^A^-t yA^ 


( 1 ) 


+ e^{W’f-a^W^+^W^ 


where A' and W' denote the spatial derivatives of cor¬ 
responding order parameters, and the coefficients fs,w> 
fes.w. and 7 are in general not independent and show 
a complicated dependence on doping and/or temperature, 
and on the mean free path. These coefficients are presented 
in Section III for the case of cuprates and iron-based pnic¬ 
tides. 

The free energy is written in the Ginzburg-Landau form 
of Eq. (1) in tbe vicinity of a critical temperature. How¬ 
ever, the critical temperatures of the transitions into 
a state with a finite W, respectively, A may be quite differ¬ 
ent. We assume that the doping level described by the pa¬ 
rameter /i is chosen in such a way (p = pc). that the critical 
temperature T^w for the non-superconducting OP W coin¬ 
cides with the critical temperature of the superconducting 
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transition Tg. This is possible because depends on fi, 
whereas Tg does not. Thus, the coefficients and hs,w de¬ 
pend on the differences rj- {\-TI Tg], d[/i^] = - /i^, and 

on impurity concentration Wimp- 

The variation of ^ with respect to A and W yields the self- 
consistency (or the Ginzburg-Lanadau) equations 

+ A[-ag + bgA^ + YW^]=0, ( 2 ) 

-iiw”+W[-a^ + b^W^ + jA^]^0, (3) 

which represent the foundation of our considerations. 

B. Homogeneous case 

Equations (2) and (3) without the spatial derivatives yield 
different uniform solutions, i.e., three different points on 
the plane of order parameters {A,W). We denote these 
points by, respectively, Ta (where A 0, W = 0), T^y (where 
A = 0, W ^0), and T(where A 0, W ^0, i.e., both OPs 
coexist), each corresponding to an extremum of the free en¬ 
ergy functional ^ (A, W). Analyzing these points we deter¬ 
mine the conditions for a particular point to correspond to 
a minimum: 

1. Ta, where A = \/Uglbg, corresponds to a minimum if 
the second derivatives of ^ with respect to A and W 
are positive, implying 

hs>0, jUg - a^bg>0. (4) 

In particular, the coefficient ag must be positive. 

2. r w, where W - \Ja^lbw, corresponds to a minimum 
if the conditions 

hw ^ 0, Yby^ ^ 0, (5) 

are fulfilled. In particular, the coefficient a^, must be 
positive. 


C. Nonhomogeneous case 


In case of heterostructures where the OPs depend on the 
coordinate x, the most interesting nontrivial solution of the 
system of equations (2) and (3) corresponds to those where 
the OP W goes to a finite value Woo = (an asymmet¬ 
ric case Woo ^ IVLoo can be considered analogously) while A 
vanishes at distances from the interface exceeding In 
other words, we consider the case when the system is at the 
point r w far away from the interface. We assume that the 
OP W(x) is suppressed near the interface, e.g., due to an en¬ 
hanced impurity scattering in the vicinity of the interface or 
diffusive scattering on the interface. The diffusive scattering 
may be caused either by interdiffusion or interface rough¬ 
nesses. As is known (see Ref. 63 and references within), the 
critical temperature T^w is suppressed by impurity scatter¬ 
ing while the critical temperature Tg of the superconducting 
transition is only weakly affected. The most essential de¬ 
pendence of the coefficients as,w and hs.w on impurity scat¬ 
tering is the one of the coefficient a^. 

Doping level near the interface also may be changed due 
to interdiffusion of atoms. Another reason for a change of 
the coefficients in the Ginzburg-Landau equations is a dif¬ 
ferent crystal symmetry at the interface. Such a mecha¬ 
nism of enhanced superconductivity at twin boundaries has 
been considered for conventional low- Tc superconductors 
in Ref. 64. A change in the energy spectrum of cuprates near 
the surface has been calculated in Ref. 65. This change can 
also lead to a modification of coefficients in the Ginzburg- 
Landau equations or even to a surface superconductivity. 
It is worth noting that the properties of surface supercon¬ 
ductivity in systems with one (superconducting) OP in mag¬ 
netic field have been studied theoretically in Refs. 66-69. 

The interface behavior of W can be modeled via the spa¬ 
tial dependence of the coefficient a^,, i.e.. 


^w(^) — 


-oq, 1x|<L, 
flw. \x\>L, 


( 8 ) 


3. Iw^, where A= y/[Ugbyi - y dw)!D and 

W = yJ(ayybg-Ya&MD with D-bgby,-Y^, corre¬ 
sponds to a minimum provided the conditions 

bg>Q, hw>0, bgbyy-Y^ >b, (6) 

are satisfied. It follows from the definition of D and 
the expressions for A and W that the conditions 

agbyy-Yaw>0, ay,bg-Y^s>0, ( 7 ) 

should be fulfilled. 

Clearly, the latter inequalities (7) are incompatible with 
those in (4) and (5). This fact is evident from topolog¬ 
ical arguments, namely, if the point Ti^a is a minimum 
of ^ (A, W), then the points Ta and T w can only correspond 
to a maximum or a saddle point of the free energy func¬ 
tional. 


where Z, is a characteristic width of the region where W is 
suppressed. The expression for aq is presented in Section III 
for a particular case. Formula for flw(^) at |x| > L implies 
that the amplitude W+oo = \/a^Jb^ is the same at x ^ +oo, 
having a lower value at the interface x = 0, see Fig. 1 (a). We 
will show that, under these circumstances, a superconduct¬ 
ing OP A arises at the interface decaying to zero as x —► oo. 

Note that in our previous publication^® we analyzed non¬ 
homogeneous solutions for the Ginzburg-Landau equa¬ 
tions (topological defects) assuming that all the coefficients 
in these equations are constant. The found solutions may 
correspond to metastable states with energies higher than 
that for an uniform solution. In the case considered here, a 
nonuniform solution for A(x) is enforced by built-in defect 
described by the a non-constant coefficient ag{x). 

In order to find the spatial dependence of W (x) and A(x), 
we assume that the superconducting OP A is small. Then, 
in the main approximation the equation for W acquires the 
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FIG. 1. (Color online.) (a) Sketch of considered system. The CDW 
(or SDW) order parameter IV is suppressed near an interface be¬ 
tween two materials in which superconducting order parameter A 
may exist alongside fV. This suppression leads to the appearance 
of interfacial superconductivity, which can he thus regarded as 
‘hidden’, (h) The case of strong suppression of W with solutions 
of A„(x) given hy hypergeometric functions. Note that Ao(jc) has 
the shape of a soliton, whereas other solutions have nodes, (c) In 
the case of a weak suppression of fV, one has a shallow ‘poten¬ 
tial’ in the ‘Schrddinger’ equation and there exists only one ‘en¬ 
ergy’ level given hy Eq. (16). 


form 


-Cw” + W[-a^{x] + b^W^]^0, (9) 

where flw(^) is given by Eq. (8). This equation can be 
solved exactly, but for simplicity we restrict ourselves to 
the simplest case of a narrow region of suppression, i.e., 
L « fw/thus obtaining the solution 

W(x) = Wootanh[)fw(|x|-i-x:o)], (10) 

where the integration constant xq obeys the equation 

sinh(2xoJCw) = 4f^jCw/aoi= f, (H) 

with jfw = f \/awl2. 

Next, we consider separately the cases of a strong [r « 1, 
cf. Fig. 1 (b)] and, respectively, weak [r » 1, cf. Fig. 1 (c)] 
suppression of W at the interface. 

Strong suppression, r«l. In this case, the product 
ZxqKv, r is small and the quantity xq in Eq. (10) can be ne¬ 
glected. Substituting W (x) in this approximation into Eq. (2) 
one obtains 


( 12 ) 


where fs = fs\/^, S-ashw-ja^, and 

g - bsb„. Equation 12 has a form of the Gross-Pitaevskii 
equation.^^’^^ Solving this equation one can determine 
the spatial dependence of the superconducting OP A. 
The calculations in this case formally coincide with those 
carried out in Ref. 70 if the interchange A « IV is done. 
For completeness we repeat the steps one has to perform 
seeking for a solution A(x). Assuming that A is small, the 
right-hand side of Eq. (12) can be neglected and we are 
left with the linearized Gross-Pitaevskii equation, i.e., the 
‘Schrddinger’ equation, which can be solved considering 
the eigenvalue problem 

^An^SnAn, (13) 

where the operator - ^cosh“^()fwX), and A„ 

and Sn are the eigenfunctions and eigenvalues of if. The 
solutions An corresponding to a discrete spectrum of Sn can 
be expressed in terms of hypergeometric functions and the 
‘energy’ levels [S < 0) being given by^^ 


’2 ,..2 


Sn - ~ 




-(l + 2n) + , 1 + 


4^ 

72^2 

Ss 


(14) 


Provided the inequality is fulfilled, there is 

only one ‘energy’ level with n = 0 and the corresponding so¬ 
lution Ao(x) has a form of a soliton. Otherwise there are 
several solutions decaying far away from the interface and 
corresponding to Sn- 

Representing the OP A close to a certain ‘energy’ level Sn 
as A(x) = c„A„(x) -I- 5An{x) with a small correction dA„(x) 
orthogonal to A„(x), one obtains the coefficients c„, 

n S — Sn 

4 =-- (15) 

g«A^(x)» 

where {(fix))) = ff^dxfix) (double angle brackets are 
used to distinguish the notation from the averaging over 
momenta directions introduced in Appendix A). 

Note an important point. The condition S <0 that de¬ 
termines the appearance of the interface superconductivity 
coincides with the condition (5) that provides the stability 
of the state with W ^0 and A = 0 in the bulk. This means 
that if a nonsuperconducting state in the bulk is character¬ 
ized by a nonzero OP W, any suppression of IV leads to the 
appearance of local superconductivity. We demonstrate this 
considering the case of a small suppression of W. 

Weak suppression, r » 1. In this case, one obtains from 
Eq. (11) XQjfw - ln\/^ and the spatial dependence of A(x) 
is determined by the ‘Schrddinger’ equation. (12) with the 
‘potential’ ^cosh“^(7CwX) ^ V(x) = 1 - tanh^ [jfwd^l + J^o)]. 
In other words, the function A(x) is determined by the 
‘Schrddinger’ equation that provides the ‘energy’ levels in a 
shallow potential well V (x). As is well known,^^ there always 
exists a single ‘energy’ level 

<fo = -//(2^s)^ (16) 


where / = {{¥ (x)>> = 4 jc“^ exp(-27CwXo). The amplitude cq is 
given again by Eq. (15) with n-Q. This means that a super¬ 
conducting condensate with a small amplitude A == Tf,Sa, 


A" +[S + ^ cosh ^ (jfwX) ] A = gA^ 
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FIG. 2. (Color online.) Sketch of temperature dependence of 
the amplitudes An corresponding to different eigenvalues Sn, see 
Eq. (14). At a given temperature T, the state with the largest A„, i.e., 
the state Aq, corresponds to a minimum of the free energy. How¬ 
ever, the transitions between different An are possible 

where Tg is the superconducting transition temperature in 
the bulk, necessarily arises at the interface as soon as the 
competing OP is arbitrarily weakly suppressed. 

Note that, in the case of a two-dimensional point de¬ 
fect instead of the interface, \So \ is an exponentially smalF^ 
quantity and, therefore, the radius of the decay of the con¬ 
densate is exponentially large. 


III. APPLICATION TO CUPRATES AND PNICTIDES 

A. Relation of coefficients to microscopic parameters 

Here, we present the expressions for the coefficients in 
the Ginzburg-Landau expansion for the case of cuprates 
and iron-based pnictides. 

As has been shown in Ref. 74, the model that has been de¬ 
veloped in detail in Refs. 75-77 for Fe-based pnictides (gen¬ 
erally, to two-band superconductors with an SDW), is ap¬ 
plicable to quasi-one-dimensional superconductors with a 
CDW, and, after certain modification, also to cuprates. 

First, we consider the region |jc| > L and assume that 
the impurity concentration in this region is small, i.e., the 
mean free path / » ^s,w- In the model of Refs. 75-77, the 
coefficients are related to the microscopic parameters of 
the model via -rj, hg - 1-05, flw = ?7(1 - Pi) - 
K = Ssm, 7 = ^ 2 m, where 77 = 1 - F/ Fg and 6[i/] ^ - nl 

with fi being a function that describes the curvature of the 
Fermi sheets of the quasi-one-dimensional superconductor 
or a deviation of the average Fermi surface from the perfect 
circle in iron-based superconductors, where partial nesting 
between the elliptical electron bands and the circular hole 
bands leads to the formation of the spin-density wave. The 
functions S 2 m, Ssm and / 3 i ,2 (see Appendix A) depend on 
the dimensionless critical curvature m = defined in 

such a way that the critical temperature Fq^ of the forma¬ 
tion of the OP W equals Fg, where Fg are the critical tem¬ 
peratures for the transition into the, correspondingly, super¬ 
conducting and CDW or SDW state in absence of the com¬ 
peting order and doping. The critical pc is determined by 


the equation (see Appendix A) 

(2p2si^(p,)>=ln(Fdw/rg), (17) 

where the critical temperature Fq^ depends, generally 
speaking, on impurity concentration which is assumed to 
be small far from the interface. 

Next, consider the region | x| < L, where the impurity scat¬ 
tering is assumed to be stronger. In this case, the temper¬ 
ature Fdw=: Fdwo[l-(47rFd„oD"^], where Fdwo is the crit¬ 
ical temperature of the transition into a state with W ^0 
in the absence of impurities and superconductivity, t -Hv 
is the momentum relaxation time for intraband scatter¬ 
ing, and the mean free path I is assumed to be larger 
than I'/Fdwo- One can easily show that, in this case, 
oq- where flw is given by the expression 

above. Provided the condition flw « (4:/rFdwoT)~^ « 1 is ful¬ 
filled, then the coefficient aq in Eq. (8) is positive (meaning 
that -flo < 0 ) and considerably esceeds flw. 

B. Temperature and doping dependence 

Considering the expression (15) for the coefficients c„, 
one can see that at small difference IS - Sn\, the OP A is also 
small. It turns to zero at certain temperatures or doping 
levels, where S[T„,fim) - <S’n(F„,p„) holds, with Sn{Tn,Hn) 
given by Eq. (14) with coefficients expressed through the 
microscopic parameters and the temperature. In Fig. 2 we 
plot the temperature dependence of the amplitudes A„ cor¬ 
responding to different eigenvalues Sn- Clearly, when the 
temperature F becomes lower than the temperature Fq de¬ 
termined by Eq. (14), the superconducting OP Ao(x) arises 
at the interface with the amplitude increasing when F is 
lowered. The temperature Fq is lower than Fg (77 > 0), but, 
under certain conditions, higher than the temperature Fb at 
which the superconducting state becomes more favorable 
in the bulk. At F < T\, a new branch Ai(F) appears, etc. 

The minimal temperature F„ (or maximal TZmax) at which 
the branch A^^j^^lF) appears, is determined by the condi¬ 
tion 277max ^ ^ 1 -t 4 ^^^^ 7 fw^ - 1. It Can be shown that at 
a given temperature F, the state with the largest A„, i.e., 
the state Aq, corresponds to a minimum of the free energy. 
However, the transitions between different A„ are possible 
analogously to transitions between an overcooled state and 
equilibrium. 

C. Interface superconducting transition temperature 

Our considerations concerned the case when the system 
is at the point F w far away from the interface, i.e., the con¬ 
ditions (5) are valid while the conditions (4) are violated. 
When applied to the case of quasi one-dimensional mate¬ 
rials with the CDW or to material like iron-based pnictides 
with an SDW, these conditions can be presented in the form 

A2T] + B2SI^^] <0, AiTJ + BiS[/I^] <0, S3m>0, (18) 
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where = ^ 3 ^ - S 2 m(l - jSi), 

A 2 = S 2 m-ssil-Pi), and iJ 2 = 105)02 if expressed via 
the microscopic parameters of the model. The coeffi¬ 
cients Ai ,2 and fJi ,2 depend on /i(/io>Aii))) and may attain 
positive or negative values (only B 2 is a positive quantity). 
If these coefficients are all positive, then this conditions 
can be fulfilled provided d[/i^] < 0. Thus, we can rewrite 
them as B 2 \d[i-fi]\ > A 2 ri and fli|d[jU^]| > Aii]. In terms of 
microscopic parameters the latter can be written as 

r]<c^Bi/Ai = l-tw, T]<c^B 2 lA 2 = l-f\, (19) 

where = |d[/i^]| and we introduced the ‘critical’ tempera¬ 
tures fw,A - Tw.aH^T^vi.s)- 

One can distinguish two from the physical point of view 
different cases; 

a) The case Ta < Tw is realized if the quantity 

D = A 1 B 2 - A 2 B 1 - S3S3m - ■sfm ^ ^ provided that Ai 
and A 2 are positive. In this case, the pure supercon¬ 
ducting state (minimum of the free energy at Ta) 
exists at temperatures T < Ta with lbs-In the 

temperature range Ti 3 ,< T < Tw, a mixed state (the 
state of coexistence) with A 0 and W ^0 given by 
expressions just before Eq. (6) takes place. At T > Tq, 
a pure CDW-state or, more generally, a W-state occurs 
with = a^ilbvi- In the interval Tw < T < Tq, 
one has a surface (or interface superconductiv¬ 
ity), where Tq is the temperature determined by 
Eq. (14) corresponding to n = 0 (the ground state). At 
To< T < Tdw, the system is nonsuperconducting with 
the OP W ^0. In Fig. 3 (a) we show schematically 
the temperature dependence of A and W and also 
the temperature range in which the local supercon¬ 
ductivity exists. At temperatures Tw,a, second-order 
phase transitions occur and the OPs A and W arise 
continuously (see expressions for A and W before 
Eq. (6) where A ~ \/Tw - T and W ~ \/T- Ta). Note 
an obvious analogy with conventional second- 
type superconductors in a magnetic field 

The quantities Ta.w are analogous to the critical 
fields Hci,c 2 [cf. Fig. 3 (c)] so that at T < Tw one 
has a purely superconducting state (full expulsion 
of the magnetic field), a mixed state in the interval 
Tw <T < T/i, (correspondingly, the Abrikosov’s vortex 
state) and a surface or interface superconductivity 
in the range Ta < F < To (correspondingly, in the 
range Hc 2 < H < Hc 3 ). At last, at T> Tq, one has a 
pure VP-state which corresponds to the normal state 
in conventional superconductors. Note that same 
analysis can be performed if, instead of temperature, 
the doping (or curvature) given by in Eq. (19) 
is considered and one obtains the corresponding 
situation as depicted in Fig. 3 (d). 

b) The case Ta > Tw is realized if the quantity 

D = A 1 B 2 - A 2 B 1 - S3S3m - < 0 provided that Ai 

and A 2 are positive. In this case, the pure supercon¬ 
ducting state (minimum of the free energy at Fa) 


exists again at temperatures T < Ta. but this min¬ 
imum is global only at T < Tipt, where the critical 
temperature Tipt for the first-order phase transition 
is determined by the equation hw) = ^{a&lbfP. 

At T> Tdw, a uniform solution for W - \/a^/bw 
arises, but it corresponds to a global minimum at 
T > Tipt. In this case, no region of coexistence exists, 
and at r = Tipt a first-order phase transition from 
the superconducting state to a state with W ^0 takes 
place with increasing temperature, see Fig. 3 (b). 
Now, if Tq as determined by Eq. (14) corresponding to 
n-0 (the ground state) falls into the region Tipt < Tg, 
which is possible if the difference Tg - Tipt is pos¬ 
itive and sufficienfiy large, then, at the interface, 
superconductivity is induced. 


D. Experiments 

The obtained appearance of superconductivity (or en¬ 
hancement of the critical transition temperature) at an in¬ 
terface between two materials in which, alongside super¬ 
conductivity, another exists and is energetically more fa¬ 
vorable, may be realized in two prominent examples of 
such systems. One of them, the cuprates, show a charge- 
density wave order alongside superconductivity^®”®® In 
these materials, an enhancement of superconducting tran¬ 
sition temperature has been found in a bilayer constructed 
of Lai.65Sro.35Cu04 and Lai.875Bao.i25Cu04.®'^^ 

Superconductivity accompanied by a spin-density wave 
is known to exist in the iron-based pnictides, where the in¬ 
terface superconductivity is proposed to be the driving ef¬ 
fect behind the almost doubling of superconducting transi¬ 
tion temperature in CaFe 2 As 2 .®®’®^ 

Unfortunately, there are no data on spatial dependence 
of the order parameter accompanying superconductivity in 
these experiments, neither the spatial dependence of the 
superconducting order parameters has been investigated in 
these experiments. Such a measurement would provide a 
test of our theory, if charge- or spin-density wave would 
have been suppressed near the interface. 

Another interesting effect also serving as a test of our 
predictions is related to the fact that there might appear a 
hysteretic behavior stemming from the presence of differ¬ 
ent ‘energy’ levels [see Eq. (14)]. This results in a sequence 
of temperatures at which A formally vanishes but the 
temperature dependence of A is determined by the highest 
temperature of them all, i.e., by Tq since the minimum of 
the free energy is deepest here. Interestingly, at a temper¬ 
ature Tn a local minimum of the free energy is present and 
adjusting A by means, e.g., of an external field, one can let 
it follow the temperature dependence of A„ after the relax¬ 
ation of external constrains, so, indeed, “there’s a lot of room 
for new combinations”.®® 
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FIG. 3. (Color online.) (a) In the case < T^, three ranges 
of temperature. For T <Ti^, the system is in a pure supercon¬ 
ducting state, whereas in the temperature range < T < Tw, a 
mixed state (the state of coexistence) with A 0 and VK 0 is real¬ 
ized. Thus, the hulk superconducting transition temperature cor¬ 
responds to Ty\f. The temperature range of the coexistence state 
is widened (enhancement of superconductivity) on appearance of 
the interface superconductivity with the highest transition tem¬ 
perature To (there may he other transition temperatures Tn)- Fi¬ 
nally, for T> To, the system is in a pure IV-state (CDW or SDW] 
up to its transition into the normal state at (h) In case 

Ta > Tw< system may he in either the pure superconducting 
or in the pure VF-state. The transition from one into another is 
of the first order at a temperature Tipt. If Tq falls into the re¬ 
gion Tipt < Ts, which is possible if the difference Tg - Tipt is pos¬ 
itive and sufficiently large, then, at the interface, superconduc¬ 
tivity is induced, (c) An analogy with conventional second-type 
superconductors in a magnetic field H is sketched. The quanti¬ 
ties and Tq are analogous to the critical fields Hd^cZ and 
respectively, i.e., up to ffcl the system is in a purely supercon¬ 
ducting state (full expulsion of the magnetic field), in a mixed 
state in the interval Hd <H< Hc 2 , and in the Abrikosov’s vortex 
state for Hc 2 < H < Hcs loosing the superconducting properties for 
H > Hc 3 . (d) If, instead of temperature, the doping (or curvature) p 
is considered, the situation is similar to the temperature depen¬ 
dence. 


IV. DISCUSSION 


We have studied a system with two competing OPs one of 
which is the superconducting OP A and another, W, can he 
the amplitude of the charge- or spin-density wave. On the 
basis of Ginzburg-Landau equations we have shown that 
if the temperature and doping are chosen in such a way 
that the state with W yiO and A = 0 is favorable, i.e., it cor¬ 
responds to a minimum of the free energy in the bulk of 
the sample, an arbitrary small suppression of W at an in¬ 
terface or a defect leads to the appearance of local super¬ 
conductivity. This mechanism of local superconductivity 


in a system with two OPs may be responsible for the in¬ 
terface superconductivity observed in many materials in¬ 
cluding cuprates and iron-based pnictides. As is firmly es¬ 
tablished, in cuprates and iron-based pnictides, CDW (or 
quadrupole charge order®^) or SDW can exist alongside with 
superconductivity. 

We found that in case of a strong suppression of W at 
some point, there are several solutions for the supercon¬ 
ducting OP A(x) which are localized on the scale of the 
coherence length These solutions are found from the 
nonlinear Schrddinger equation (or Gross-Pitaevskii equa¬ 
tion) and correspond to different ‘energy’ levels. Each solu¬ 
tion arises at certain values of temperatures Tn (or doping 
level and has a different form changing from a soliton- 
like one to an oscillatory function decaying at infinity. How¬ 
ever, only the soliton-like solution Ao(x) corresponds to a 
minimum of the free energy. Other solutions A„(x) (with 
n^O] with nodes have higher energies. They correspond to 
metastable states. If the corresponding potential Y (x) in the 
Schrddinger equation is not deep enough, there is only one 
‘energy’ level and only one soliton-like solution for A(x). 

In the case of an asymmetric potential Y (x), the localized 
solution for A(x) exists provided that the potential well is 
deep enough.^^ If for some reasons the OP W is suppressed 
at the surface, the solutions for A(x) have the same form as 
in case of the symmetric Y (x) and surface superconductiv¬ 
ity arises in the sample. Note an important point. The lo¬ 
cal superconductivity may arise in one- or two-dimensional 
cases (a flat interface or a point defect) if the suppression 
of the OP W is small. This means that even if the minimal 
value of W(x) corresponds to a minimum of the free energy 
in an uniform case, in a nonuniform case local supercon¬ 
ductivity would arise at the interface or at the surface. 

The developed theory is able to explain the emergent or 
enhanced interface superconductivity in some cuprate or 
iron-based pnictide heterostructures. The predictions can 
be tested in further experiments. 

Remarkable, the used approach based on the consider¬ 
ation of Ginzburg-Landau equations for two coupled or¬ 
der parameters does not depend on the nature of order de¬ 
scribed by those. The obtained results are also applicable to 
description of an arbitrary ‘hidden’ order evolving at the in¬ 
terface, e.g., the ferromagnetism induced at the interface of 
an oxide heterostructure observed recently.^^ 

Note that localized superconductivity may arise also in a 
homogeneous sample if a nonuniform W (x) state (for ex¬ 
ample, stripes) is energetically favorable due to an inter¬ 
nal mechanism (e.g., the Larkin-Ovchinnikov-Fulde-Ferrel 
mechanism).Analysis of this state with two competing 
order parameters deserves a separate consideration.®® 
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Appendix A: Doping dependence of the coefficients in 
Ginzburg-Landau equations 

In the notation of Refs. 75-77, the Ginzburg-Landau 
equations have the form 

+ A[wh2m + - mTJT)] = 0 , 

(Al) 

-e^V^W+W[{2^J.hlm) + Wh3m + ^^S2m-mTa^lT)]^0, 

(A2) 

where ^s,w are the coherence lengths (at low temperatures) 
for A and W, respectively, and Ps.dw are, respectively, the 
critical temperatures for the transition into the pure super¬ 
conducting state or into a state with a CDW or an SDW 
only. In other words, Pdw is the critical temperature for 
the transition into the charge-ordered state in absence of A 
and [I, while Ts is the superconducting transition temper¬ 
ature in absence of W. The angle brackets mean the an¬ 
gle averaging (in Fe-based pnictides) or integration along 
the sheets of the Fermi surfaces in quasi-one-dimensional 
superconductors. The functions sim, S 2 m> etc. are func¬ 
tions of the normalized curvature m- ^l{nTs), where 
/i = fio + ^ipCos[{p^ + is a curvature in quasi-one- 

dimensional superconductors with a doping-dependent 
value of Pq. It is assumed that the Fermi surface of these su¬ 
perconductors consists of two slightly curved sheets which 
are perpendicular to the x axis^^. In the case of Fe-based 
pnictides, p = po + pq, cos(2^) is a quantity that describes an 
elliptic {pq, 0) and circular [pq, = 0) Fermi surfaces of elec¬ 
tron and hole bands. All quantities—A, W and p —are 

measured in units of ttTs- The expressions for the coeffi¬ 
cients in the G-L expansion with account for impurity scat¬ 
tering have been calculated in Ref 86. 

Replacing the derivative V ^ V - i2:;rA/<l>o. one can use 
Eqs. (Al) and (A2) to describe vortices in superconductors 
with a CDW,®^ where <I>o is the magnetic flux quantum. 

As it is seen from Eq. (A2), the critical temperature 
depends on doping, i.e., on the parameter p. We choose 
this parameter p = pc in such a way that = fs- This 

means that at T = T^, the quantities A = W = 0, and, thus, pc 
obeys the equation 

(2pcSim(A«c)> =ln(rdw/rs), (A3) 

where pc is a function of two parameters, 
i.e., pc = Pc{po,Pq,). 

Then, we expand the function s\m{p,T) in the de¬ 
viations 5[p^]-p^-p^ and ST-T^-T, thus obtaining 
Simip, T) = Simipc, Ts) + pi6T + {P 2 S[p^]), and use Eq. (A3) 
to obtain equations in a general standard form (assuming 
that all the functions depend only on one coordinate x), 

+ A[-as + bsA^ + rW^]^0, (A4) 

-^iw”+W[-a^ + by,W^ + jA^]^0, (A5) 


with A' and W' as well as A" and W" denoting the first and 
second derivatives with respect to x, respectively. These 
equations determine extrema of the free energy functional 
[cf Eq. (1)] 


W' 


/ dx 


^2(A')2-asA2-t yA^-tyA^W^-t (A6) 




with respect to A and W, and the corresponding coefficients 
of the G-L expansion are related to variables in Eqs. (Al) 
and(A2) via as -T],bs- S 3 =^ 1.05, av,-r][l- - (P 2 S[p^]), 

bw = -Ssmi 7 = S 2 m, where 77 = 1 - T/Tg. The expressions for 
the coefficients in terms of the microscopic parameters of 
the model for cuprates and iron-based pnictides are given 
as follows: 

00 

S3=X(2« + l)■^ (A7) 

n=0 

00 

Y^i2n+ir^[(,2n+lft^ + m^]~^, (A8) 

n=0 

S2m = Y. ([(2n+ 1)^- m^](2n+ 1)"^[(27?+ ifm^]~ ), 

n=0 

(A9) 

00 

S 3 m= E {i2n+l][{2n+lf-3m^][{2n+lf + m^]~^), 

n=0 

(AlO) 

(4m^(277+l)[(2/7+l)^ + m^] ), (All) 

n=0 

00 

782 = E 2{2n+iy'^[i2n+lf + m^]~^, (A12) 

n=0 

where t- T/Tg, the angle brackets (...) denote the angle av¬ 
eraging (in iron-based pnictides) or integration along the 
sheets of the Fermi surfaces (in quasi-one-dimensional su¬ 
perconductors or cuprates). 


Appendix B: Details on solution of the Gross-Pitaevskii equation 


Here we sketch the solution of the Gross-Pitaevskii equa¬ 
tion for A. In zero-order approximation we obtain for Aq 
from Eq. (12) 

|^^Ag-i-Ao[<S‘-i-'^cosh“^(XwX)] =0. (Bl) 

This equation is integrable and its solutions correspond¬ 
ing to a discrete spectrum of Sn are expressed in terms of 
hypergeometric functions^^. In our notations, the ‘energy’ 
levels of discrete spectrum are given by^^ 


Sn — 


r 2„2 


-(l + 2n) + , 1 + 


4^ 

2^2 




(B2) 


and their maximal number ttmax is determined by 
2 ainax — ^1 + 4‘^^g — 1 . 
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We expand the correction SA to the zero-order solu¬ 
tion Ao in terms of the normalized eigenfunctions A„ = i//„ 
of the operator ^ ~ ^cosh“^(jCwX). These func¬ 

tions obey the equation 

■^y^n — ^nV^n- (h3) 


Solutions of Eq. (12) can be written explicitly if the quan¬ 
tity S = S(ri,6lfi^]] is close to a certain ‘energy’ level Sn, say 
to Sn, such that S - Sn = S{riN,Slfi^^]) (if considering the 
model for cuprates or iron-based pnictides, the ‘tempera¬ 
ture’ T] or doping 6 [p^] should be chosen properly). We write 
Eq. (12) in the form 


^A^SnA + R{A), (B4) 


with R- gA^ + {S -Sn)A and represent A as 
A(x) = CNy/N(x) + SAnM, where SAn(x] = CN.nVnix), 
and the summation runs over all n except the term n- N. 
We substitute this A(x) into Eq. (B4) and multiply this equa¬ 
tion first by yfN and then by y/n with N, then integrating 
the obtained result each time over x. Thus, taking into 
account the orthogonality of different eigenfunctions, we 
find the coefficients c„ 


% - 


S — Sn 


8((Vn'>'> 


CN,n — S^N 


3 ((V%Vn)) 


Eti-En 


with N, 


(B5) 

(B6) 


where ((/(x))> = /^dx/(x). Obviously, in Eq. (B6), yjn 
and yiN have to have same parity (both even or both odd). 
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